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Abstract: We update the allowed parameter space of the simple chaotic inflationary 
model with quartic potential and light inflaton [1] taking into account recent results from 
cosmology (CMB observations from SPT, ACT and WMAP) and from particle physics 
(LHC hints of the SM Higgs boson). The non- minimal (yet small) coupling to gravity of 
the inflaton becomes essential to fit the observational data. The inflaton has mass above 
300 MeV and can be searched for at S-factories in S-meson two-body decays to kaon 
and inflaton. The inflaton lifetime depends on the model parameters, resulting in various 
inflaton signatures: either a missing energy, or a displaced vertex from the S-meson decay 
position, or a resonance in the Dalitz plot of a three particle decay. We also discuss the 
implementation of the inflaton model to the z^MSM, where the inflaton can be responsible 
for production of the dark matter sterile neutrino in the early Universe. 
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1. Introduction 

Elegant solutions to the major problems of the Hot Big Bang theory are provided by 
a preliminary inflationary stage of the Universe expansion [2, 3, 4, 5, 6]. The simplest 
realization is made using a scalar degree of freedom, inflaton, which evolves very slowly at 
the inflationary stage providing nearly constant energy density and hence the accelerated 
expansion. At later times the inflaton starts to oscillate and the Universe exits inflation 
and enters a stage when rapidly evolving homogeneous inflaton dominates the total energy 
density. The energy has to be transferred to the Standard Model (SM) particles thus 
reheating the Universe. To make the transfer an inflaton-to-SM coupling is needed. This 
coupling takes care of the SM particle production and the onset of Hot stage of the Universe 
expansion. 

Let us analyse the consequences of adding explicitly just one scalar inflaton field to 
the SM. In a renormalizable model of the inflaton sector the coupling to the SM should 
be renormalizable as well. Then for the simplest inflationary models with only one scalar 
field, inflaton, involved in the inflationary dynamics, the only suitable coupling to the SM 
is its coupling to the SM Higgs field via the scalar potential. Then reheating happens via 
the Higgs boson production by oscillating inflaton field. A variant of these models has been 
considered in Refs. [7, 1]. It is distinguished by the absence of dimensionful parameters in 
the sector of SM fields. The vacuum expectation value of the SM Higgs field v appears as 
a nonzero vacuum expectation value of the inflaton field: the quartic coupling between the 
two field is responsible for both reheating in the post-infiationary Universe and electroweak 
phase transition later at the Hot stage. The only dimensionful parameter is in the infiaton 
sector, and it is of the order of electroweak scale. All model parameters are therefore stable 
with respect to the finite quantum corrections (or logarithmic renormalization running). 
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The appealing phenomenological feature of the inflationary model described above is 
the new light particle, inflaton, with the light mass in GeV region. Cosmology limits 
the mass range of the light inflaton from above and from below: lighter inflaton causes 
very late reheating, while quantum corrections from the inflaton-to-Higgs quartic coupling 
spoil flatness of the inflaton potential for heavy inflaton masses. The latter bound on the 
inflaton-to-SM coupling constrains the inflaton mass to be below 2 GeV [1] . 

In the model the quartic coupling induces the Higgs-inflaton mixing, so the inflaton 
mass is proportional to the mass of the Higgs boson and square root of the quartic coupling. 
The inflaton mass range presented above refers to the Higgs mass of 100-200 GeV. Recent 
results from LHC [8, 9] can be interpreted as evidence for the SM 125 GeV Higgs boson, 
which allows to resolve the ambiguity due to the Higgs mass value: now the only free 
parameter remaining in the model is the Higgs-inflaton quartic coupling, or equivalently 
the inflaton mass. Thus the estimates given in Ref. [1] can be appropriately refined. It 
is important for the phenomenology of the light inflaton, studied there. Indeed, provided 
mixing with the SM Higgs boson, inflaton may be produced in scatterings and decays of 
the SM particles and subsequently decays into the SM particles. The latter branching 
rations are the same as those for a (hypothetical) SM Higgs boson of mass equal to those 
of the inflaton. Thus, the inflaton mass alone determines the inflaton production rate in 
the laboratory experiment and its lifetime, which was estimated as 10~^-10~^'' s. The most 
promising place to search for the inflaton was suggested to be in B-meson decays [1] . 

Another issue has risen due to the recent analyses of cosmological data, mainly obser- 
vations of the Cosmic Microwave Background (CMB) anisotropy by SPT [10], ACT [11], 
and the common global fit together with the analysis of CMB data collected for 9 years 
by WMAP [12]. The most important, relevant for our model result is lowering an upper 
limit on the tensor-to-scalar perturbations ratio, which excludes the simple chaotic inflation 
with quartic potential discussed above at more than 95% confidence level. It was already 
pointed in Ref. [1] that with a small non-minimal coupling to gravity the situation may 
be cured. It was found [1] that, at rather small value of non-minimal coupling constant 
< 10~^, this also leads to the increase of the infiaton self-coupling constant. The latter is 
determined from the amplitude of the scalar perturbation power spectrum obtained mostly 
from analyses of CMB anisotropy. This change was accounted for in the estimates pre- 
sented in Ref. [1]. However, with new results from CMB anisotropy [10, 12, 11] a somewhat 
larger non-minimal coupling is needed to make the inflaton model viable, compared to the 
value adopted in Ref. [1]. Hence that analysis has to be modifled appropriately. 

In the present paper we update the analysis of Ref. [1] to account for both recent 
measurement of the Higgs boson mass of about 126 GeV [8, 9] and recent new upper limit 
on the tensor-to-scalar ratio r < 0.13 (at 95% C.L. [12]). 

The paper is organized as follows. In Sec. 2 we present the model Lagrangian and 
discuss the particle spectrum in the scalar sector. In Section 3, which is devoted to in- 
flationary dynamics of the model, we give predictions for cosmological parameters to be 
tested in future cosmological experiments. Particle physics phenomenology is discussed in 
Sec. 4, where we reflne the estimates of the mass range of the light inflaton and discuss the 
strategy of searches for light inflaton in particle physics experiments, such as LHCb and 
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Belle-II. Section 5 contains conclusions and further discussion on the model extension to 
solve the major phenomenological problems of the SM: neutrino oscillations, dark matter 
phenomenon and baryon asymmetry of the Universe; here an example is z^MSM [13, 14]. 



2. Model: inflaton with /3X^/4 potential, non-minimal coupling to gravity 
and mixing with the SM Higgs boson 

The action of the light inflaton model is [1] 



SxSM 



■■ J d^x (£sM + ^XN + >Ccxt + >C 



grav; ) 



CxN = Id^Xaf^X + ^mlX^ - - ^ {h^H - jX^f , (2.1) 

■i-grav — 2 \ ) 

where jCsm is the SM Lagrangian without the Higgs field potential, and Cext stands for an 
extension of the SM capable of explaining all the major phenomenological puzzles: neutrino 
oscillations, dark matter phenomena, baryon asymmetry of the Universe (see discussion in 
Sec. 5). As far as we assume a,/3 <^ A, inflation proceeds along the flat direction of the 
scalar potential, where the Higgs and inflaton contributions in the last term of (2.1) cancel. 
In the inflaton sector of the scalar potential the sign of the quadratic term is chosen to lead 
to the nonzero vacuum expectation value for the inflaton field at late stages of the Universe 
evolution (after reheating). Then the last term in (2.1) gives rise to spontaneous breaking 
of electroweak (EW) symmetry, and the SM Higgs field gains non-zero vacuum expectation 
value as well. The two scalar excitations above the vacuum, h and Xi correspond to the 
SM Higgs boson and inflaton particles, but the mass basis is slightly rotated as compared 
to H and X, with the small mixing angle 9. Four parameters of the model, mx, /3, A, and 
a, determine the Higgs field vacuum expectation value v ~ 246 GeV, the Higgs boson mass 
m/i ~ 126 GeV, and the inflaton mass 




^h\hr- = (2.3) 



Thus, at a given value of /3, the only free parameter in the scalar sector is the mixing 
coupling a or the inflaton mass m^^, which governs the inflaton effective coupling to all the 
SM flelds via mixing with the Higgs boson characterized by the squared mixing angle 

0^='-^ = '-^. (2.4) 
A 

For zero or very small non-minimal coupling constant (^ < 10""^) the inflationary dynamics 
is fully determined by the parameter /3, which is then flxed from the amplitude of primordial 
density perturbations. This amplitude is measured from observations of CMB anisotropy, 
and the resulting value, used in Ref. [1], is /3 = (1 — 2) x 10~^^. In this case the only 
free parameter in the model is the inflaton mass (uniquely related to a for given Higgs 
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mass and /3, see (2.4)). Further, stability of the inflaton potential with respect to quantum 
corrections from the SM sector due to inflaton mixing to the Higgs field places an upper limit 
on a, while lower limit comes from the requirement of efficient reheating to the temperature 
exceeding 100 GeV (to allow for some baryogenesis mechanism). Consequently, using the 
bounds for a, the inflaton mass was conflned in [1] between 30MeV and 2 GeV. In Sec. 3 
we extend those estimates for larger (^ > 10~^) values of non-minimal coupling constant 
favored from the combined analysis of present cosmological data [12]. 

Let us note that generally one more mass parameter in the SM Higgs sector may be 
expected (corresponding to the quadratic term for the Higgs field) and one more mass 
parameter in the infiaton sector (the one in front of the cubic term). The latter is even 
welcome, as it helps to get rid of domain walls emerging in the Universe at the moment 
when infiaton gains non-zero vacuum expectation value and Z2 symmetry of (2.1) becomes 
spontaneously broken. In case the new dimensionful parameters are of the order of EW scale 
or below, they do not change the scale of infiaton mass, but relaxe the relations between 
inflation- related and low-energy physics parameters: the value of inflaton mass would not 
uniquely determine the inflaton couplings to other fields in that case. In what follows we 
assume those parameters are much smaller than the EW scale and hence negligible for our 
analysis. 



3. Inflation and reheating: prediction of cosmological parameters 

The infiation in the model happens along the direction H^H = jX'^, so that the potential 
term in (2.1), which has a not quite small coupling constant A, vanishes.^ The mass term 
in (2.1) is negligible because mx is significantly below the infiationary scale. Thus, the 
only relevant terms during the inflation are the quartic potential term and non-minimal 
coupling with gravity. This brings us to the situation thoroughly analyzed in literature 
(see, e.g. [15, 16]). Below in Sec. 3.1 we repeat the analysis to arrive at already known 
results and obtain some new analytic approximations valid in the interesting part of the 
parameter space. 

3.1 Slow^ roll parameters and spectral indices 

To apply the standard slow roll formalism (see e.g. [17, 18]) we perform a conformal trans- 
formation of the metric 

9^^u ^ g^.u = ^^g^u, n^ = l + iX'^/Ml , (3.1) 

In the new variables g^u (Einstein frame, as opposed to the original Jordan frame) gravity 
couples minimally to all matter fields including inflaton. Then the inflaton potential gets 
modifled and reads 



^This is true for /3 ^ 1, which is the case we are interested in this article. 
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After the conformal transformation the inflaton kinetic term gets modified too. In the 
Einstein frame the canonically normahzed "inflaton field" X is obtained from the solution 
of 

dx = y — — • ^^-^^ 

Then the potential term for the field X is (3.2) where the solution X = X (X) of eq. (3.3) 
is substituted. 

For the canonically normalized field X we can use the standard formulas for the slow 
roll parameters during infiation, number of e-folding and spectral indices. It has been 
proven for our model that prediction for cosmological parameters does not depend on the 
choice of frame [16] . In the Einstein frame we have a usual large field chaotic inflation with 
the potential (3.2), (3.3). In this case the amplitudes and spectral indices of scalar and 
tensor perturbations, prediction of infiation, are generally determined by the value of the 
scalar potential and values of the slow roll parameters. The slow roll parameters are [15] 

MlfdU/dXy^ 8 

2 \ U J (1 + e(l + 6C)X2/M2 )X2/M2 ' ^ ■ ^ 

_ , £u/dx^ _ 4(3 + e(i + mxyMi - 2e{i + GQXVmI,) 

^ U (l + e(l + 6e)X2/M2)2x2/M2 • ^•'•^^ 

At the end of inflation e = 1 and then from (3.4) the field value X^ is given by [15] 

X 2 _ v^l92^2 + 32^ + 1 - 1 



M2 2^(1 + 60 



(3.6) 



During inflation the inflaton quantum fluctuations become a source of matter perturbations. 
We are interested in the moment during inflation when perturbations of WMAP pivot scale 
exit the horizon due to stretching in almost exponentially expanding Universe. To estimate 
the inflaton fleld value X^ at the moment, corresponding to e- foldings before the end 
of inflation, we should invert the relation [15] 

Jx^ MldU/dX\dx) """^ Jx^ 4(l + eX2/Af2)M2 



i + i") [Xl/Ml - Xl/Ml) - In 



6y' ^^' ^ '^' ^^ i + ex2/M' 



(3.7) 



Solving this transcendental equation for Xf^ analytically is impossible. To obtain accurate 
predictions we solve it numerically and present the results in Fig. 1; below we describe also 
the useful analytical approximations. 

Now we have all formulas to match parameters of scalar and tensor perturbation spec- 
tra probed in cosmological experiments. We start with the perturbation amplitudes. One 
adopts the WMAP normalization of matter perturbation power spectrum [12] to fix the 
ratio 

U/e = 2Att^A^mI, ~ (0.0276 x Mp)^, (3.8) 
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that in our model (3.2), (3.4) is 



u ^ Pirn + ml/Ml + i)x%/M% 

and allows to determine the value of /? at given ^ and as 

32(1 + /m2)2 
^ '^{i{Qi + l)X%/Ml + l)X%/Ml,' ^ > 

For numerically solved (3.7) at relevant = 60 (see below Sec. 3.2), we present /3 as 
function of in Fig. 1, lower right panel. 

Let us proceed with tensor-to-scalar ratio r and spectral indices of scalar and tensor 
perturbations (n^ — 1) and ut, respectively. To the leading order in slow roll parameters 
(3.4), (3.5) the general relations are linear (see e.g. [17, 18]): 

r = 16e, n, -l = 2r/-6e, nr = -2e. (3.11) 

They must be evaluated at the horizon crossing, when X = Xjy, see eq. (3.7). 

Two opposite cases are easy to investigate: models with small and large values of ^. 
Indeed, in the limit ^ — )• we are back to the chaotic inflation with quartic scalar potential 
[19]. Then from eqs. (3.6), (3.7) Xj^/Mj, 8{N + 1) and further from eqs. (3.4), (3.5): 

13 , , 



N + V ' 2{N + l)' 
Plugging these results into eq. (3.11) one arrives at 



16 3 2 , ^ 

ns — ^ = -, TiT = —— -. (3.13) 



N + 1' N + 1' N + 1 

Eq. (3.10) gives for the quartic self-coupling 

^TT^A^ (3.14) 
^ 2(A^ + 1)3' ^ ^ 

which equals 1.5 x 10^^^ for the relevant value N=60 (see Sec. 3.2). It was adopted in 
Ref. [1] as a reference value for all numerical estimates. 

In the opposite limit ^ — t- oo one obtains from eqs. (3. 6), (3. 7) at large N (i.e. neglecting 
logarithmic correction in (3.7)): S^X'^/Mp — )• 4A^/3 and further from eqs. (3. 4), (3. 5): 

^ 1 

^ = -1^- (3-15) 



47V2' 

Plugging these results into eq. (3.11) one arrives at [20, 16] 



12 2 3 , , 



Self-coupling depends on as follows from eq. (3.10): 



72^2 A2 „ 

/5 = ^^C^ (3.17) 
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that is 13 ^ 4.5 x 10-^°^^ foj. relevant = 60. Thus even /3 ~ 1 is allowed at sufficiently 
large ^, and hence the SM Higgs boson non-minimally coupled to gravity may play a role of 
inflaton [21]. Note that in the limit of large ^ formulas (3.16) coincide with those obtained 
within i?^-inflation [2, 3, 22]; however, numerical predictions generically differ because of 
different reheating temperature resulting in different values of N in different models, see 
discussion in Sec. 3.2. In particular, the numerical predictions of i?^-inflation [3, 22, 23] 
and inflation driven by the Higgs boson [24] differ [25, 26]. 

To find a reasonably good analytical approximation for the intermediate case of small 
but finite value of ^ we utilize the following relation 

(1 + 60 ^|^/M|, = 8 (iV + 1) , (3.18) 

that follows from eqs. (3.6), (3.7) when linear in ^ corrections are included. Quite re- 
markably, though formally obtained at small ^, eq. (3.18) interpolates smoothly between 
—7- 0" and —7- oo, large iV" regimes considered above, with only \nN/N ~ 5% relative 
error in in the large ^ limit. Putting (3.18) into (3.4), (3.5), and (3.11) one obtains 
formulas for tensor-to-scalar ratio r and spectral indices, interpolating between (3.13) and 
(3.16), cf. [27]: 

16(1 + 60 



(iV + l)(l + 8(iV + l)0' 



(3.19) 



^ _ 3(1 + 60 + 8(iV + 1)(5 + 24g)^ + 128(jV + l)^g^ 

(iv + i)(i + 8(iv + i)e)2 ' ^'^■^^> 

2(1 + 60 

"^-"(iV + l)(l + 8(iV + l)0- ^^■^^> 
Similarly one obtains from (3.18) and (3.10) the formula for inflaton quartic coupling, 

^ Stt^ A| (1 + 60(1 + 6g + 8(jV + 1)0 
^ 2 (l + 8(iV + l)0(iV + l)=^ ' ^ ' 

interpolating between (3.14) and (3.17). 

The approximations (3.19)-(3.22) are functions of ^ and N. In our model with natural 
reheating through the Higgs boson production we obtain = 60, see Sec. 3.2. In Fig. 1, 
we present n^ — 1, r and /3 as numerical solutions. The approximate formulas (3.20), (3.19), 
and (3.22), have relative errors of not more than 8%, 4%, and 8%, respectively, in the entire 
interval of (and are much more precise at low 0- From Fig. 1 one concludes that the 
cosmological parameters effectively reach asymptotes of zero and infinite C ^ ^ 10~^ and 
^ > 0.1, respectively. 

Note finally, that the analysis in this section can be spoiled by the radiative corrections 
from the interaction with the Higgs boson, which gives the contribution of the order to 
the inflaton quartic coupling. We will somewhat arbitrary use the bound of < 0.1/3 for 
the radiative corrections to the inflaton. 

3.2 Horizon crossing and reheating 

Reheating in the model happens via the SM Higgs boson production through inflaton-to- 
Higgs coupling in (2.1). Usually the inflationary predictions are plotted depending on the 
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0.93 0.94 0.95 0.96 0.97 0.98 0.99 1 10"^ lO"'* 10"^ 10"^ 10"^ 10° 

Figure 1: The cosmological parameters for < ^ < 1; number of e-foldings is set to = 60 (see 
discussion in Sec. 3.2). It is clearly seen that already at ^ = 0.1 the model reaches the asymptote of 
"large" non- minimal coupling. The shaded regions on the (r, n^) plot are outlined by la- and 2(7- 
bounds from WMAP9 data (green) and WMAP9+BA0 data (red) [12]. The interval ^ < 2 x 10"^ 
is disfavored from cosmology. 

number of e-foldings N that happened after inflaton perturbations of a specific observable 
scale exit the horizon. This number is commonly estimated as = 50-60. Let us stress, 
that in the model of this article the dynamics of the Universe expansion is known for the 
whole time after inflation, including reheating. Specifically, if the non-minimal coupling 
is not too large, post-inflationary dynamics is dominated by the quartic scalar potential. 
This means that the Universe expands as at the radiation dominated stage. The actual 
reheating moment, though happening quite late [27], leads only to redistribution of energy 
between different relativistic degrees of freedom, and does not influence the rate of the 
Universe expansion. 

Let us derive the relation between the number of e-foldings N and model parameters. 
We are interested in the moment when the perturbation modes of conformal momentum 
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k, corresponding to the WMAP pivot scale k/ao = 0.002 Mpc~^ [12], exit the horizon at 
inflationary stage. Hereafter oq is the present scale factor, "e" and asterisk in subscript 
refer to the values at the end of inflation and at the horizon crossing of the modes deflned 
above; so iV = log(ae/a*) by deflnition. At horizon crossing the Hubble parameter Ti obeys 
the equation 

Then 

iV = iog^ + log^+log^, (3.23) 

where "r" in subscript refers to the values at reheating. At the hot stage entropy in 
comoving volume is conserved, hence 

rr 1/3 

^_n9^ f3 241 

where for the effective number of the SM degrees of freedom at present time and at reheating 
one has go = 43/11 and gr = 106.75, respectively. Since between inflation and reheating 
the Universe in our model expands as it would be at radiation domination, we have 

rr 1/4 1/2 



ar ul'^ 301/4 



(3.25) 



Plugging eqs. (3.24), (3.25) into (3.23) and substituting the Hubble parameter from the 
Friedman equation we obtain 

7^ 1/3 1 2 rrl/2 

N = log ^ + J log + log (3.26) 

k/ao 4 270 MpUl'^ 

where Tq is CMB temperature at present. Note that dependence of on the reheating 
temperature enters only through the number of degrees of freedom at reheating g^-. It is 
due to the fact that at hot stage the expansion is not isothermal but an adiabatic process, 
so the cosmological observables depends on the effective numbers of degrees of freedom in 
plasma at the onset of the hot stage gr and at present g^. 

Numerically, for WMAP pivot scale k/ao = 0.002 Mpc"i [12], one gets 

r7l/4/),T rrl/4 1/4 

A = 64.3 + log — + log . (3.27) 

1 /4 

The flrst logarithmic term in (3.27) follows from eq. (3.8), log ^* J,^^ = —3.6. Substituting 
eqs. (3.4), (3.2), (3.6) into the last logarithmic term in (3.27), one arrives at 

1 (1 + 16C+ \/l + 32(l + 6nn Aat/Mp 

AT = 60.7+ -log ^ — I — — • (3-28) 

2 sV2^l + {l + Qi)iXl/Ml{l + iXl/Ml) 

We use this formula for the numerical estimates of the number of e-foldings in a model 
with non-zero flnit ^. One gets A^ = 61.5-58.5 for ^ changing from zero to inflnity. 
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Let us analyze briefly the reheating itself. After the inflation the action for the scalar 
sector in the Einstein frame is 




/M|. d^H^d^H n'^ + eexyu^.d^xd^x v{h,x) \ 
\ 2 2 J- 



In the reheating regime X < for ,f < 1 we can neglect the coefficient in front of the 
kinetic term for X: the Universe is dominated by the field X oscillating due to the quartic 
potential. For zero ^ the reheating due to the potential mixing term aX'^H'^H was analyzed 
in detail in [27], and leads to the bound 

a > 0.7 X 10"^\ 

For ^ 7^ additional contributions to the scattering process appear suppressed^ by at least 
Mp/,^. It is easy to estimate that this contribution is negligible compared to the quartic 
mixing term and all the bounds on a in [27] remain valid. 

4. Phenomenology: searches in particle physics experiments 

We see that models with very small non-minimal coupling ^ < 2 x 10"'^ are disfavored from 
cosmology (upper panels in Fig. 1) and hence the relevant for low energy phenomenology 
quartic self-coupling /? should obey /? > 3 x 10~^^ (lower right panel in Fig. 1). Accounting 
for ^-dependence of /3 in eq. (2.3), where the SM Higgs boson mass is m/^ = 126 GeV, and 
recalling that parameter a is bounded from above and below by constraints on quantum 
corrections and reheating temperature, respectively, as explained in Sec. 2, we end up with 
cosmologically motivated interval of the inflaton mass m^, confined between green and blue 
shaded regions of Fig. 2 at ^ > 2 x 10^^. 

The upper limit on at a given ^ comes from requirement of sufficiently energetic 
reheating, that develops in the same way as in the case of = 0, studied in Ref. [27]. 
It implies a limit a > 0.7 x 10~^^, which is substituted to (2.3) and with account of ^- 
dependence of f3 is plotted in Fig. 2. The lower limit on niy. for the case of = has been 
obtained in [1] from smallness of quantum corrections to the infiaton potential, guaranteed 
at 

< 0.1 X ^. 

While ^ < 1, that estimate remains valid. We plug it into eq. (2.3) and recall ^-dependence 
of P to depict the upper limit on inflaton mass in Fig. 2. 

Another set of bounds on the model parameters comes from particle physics experi- 
ments. Thanks to the inflaton-Higgs mixing (2.4), inflaton can be produced in the same 
channels as the SM Higgs boson would be of the same mass m^. The same mixing is 
responsible for the inflaton decay to the light SM particles, and its decay pattern is exactly 

^Let us note here the significant difTerence from the situation of ^ 3> 1, wliicli is similar to tlie case of 
inflation [23, 26]. For ^ 3> 1 the kinetic term us significantly modified, and the field redefinition for X 

is needed during reheating, leading to a much weaker suppression of the interactions with the Higgs boson 

by first power of Mp . 
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Figure 2: The bounds on the inflaton in the planes (to^, ^) and {m^, 9"^). Cosmological constraints 
are shaded: the green region leads to insufhcient reheating, the blue region gives large radiative 
corrections to inflaton potential; the region ^ < 2 x 10~^ is disfavored from analysis of CMB data 
presented in Sec. 3. Other constraints are from direct searches in accelerator experiments: the 
brown region is excluded by the CHARM experiment. On the upper two plots other lines are the 
bounds from B meson decays, with the excluded areas to the right of the lines. On the lower plots 
the lines correspond to inflaton lifetime. Dashed lines on the right plots are isolines of constant ^. 



the same as that of the normal SM Higgs boson, if its mass would be equal to m^. Inflaton 
decay branching ratios are calculated in Ref. [1] and here we reproduce them in Fig. 3. 
Note that in the inflaton mass range 1-2 GeV the estimates are rather vague because of 
QCD-uncertainties with hadronization. The inflaton lifetime depends on (3 and hence on 
^ given the relation in Fig. 1 (right lower panel). We have recalculated inflaton lifetime by 
replacing /3 — ?■ /3(^) in formulas of Ref. [1], the result is plotted in Fig. 3 (right panel). 

Inflaton-to-Higgs mixing opens a room for direct searches for the light inflaton, which 
has been thoroughly studied in Ref. [1]. Remarkably, only and (3 enter inflaton-Higgs 
mixing (2.4), hence no dependence on the Higgs boson mass in the inflaton production rate 
or its decay branching ratios. Thus laboratory experiments allow for probing directly the 
inflaton quartic coupling responsible for the inflationary stage in the early Universe. In [1] 
we have found as most relevant the analysis of the CHARM bound. In case of non-zero ^ 
it goes along the same lines as in Sec. 6 of Ref. [1], with the only substitution of /3 — )• /3(^) 
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Figure 3: Inflaton decay branching ratios {left plot) and inflaton lifetime {right plot); theoretical 
predictions for ~ 1 — 2 GeV (thin dashed lines on the left plot and dotted lines on the right 
plot) suffer from significant QCD-uncertainties. 

in the inflaton decay length and meson decay branching ratios into inflaton. The resulting 
bound on is read from Fig. 5 of Ref. [1] and thus we obtain the lower limit on the 
inflaton mass at given ^ depicted in Fig. 2. One observes in Fig. 2, that with increasing 
^ at first the lower masses become allowed, but quite fast the radiative corrections bound 
for the inflation starts to dominate. Moreover, the whole inflaton window moves to higher 
masses (that was in fact known, see [27]). 

The most promising processes in searches for the inflaton of < 5 GeV are 6-quark 
decays. In Ref. [1] we found for the i?-meson decay rate 

MB - .X, . 10- X (l - <J (,^] (?^y ,4.1, 

~ 4.8 X 10"*^ X ^1 - 

Here Xg refers to a strange hadron, which most probably turns out to be a A'-meson, 
because the other decay product, x, is a scalar. For a not so small value of ^ the inflaton 
is short lived (see Fig. 3) and decays inside the detector into SM particles, 

^ ^ ^+^- , vr+TT- , TrVO , K+K- , . . . , 

thus contributing to the corresponding three-body decay modes of i?-meson. Conserva- 
tively, we require that inflaton contributions, estimated with eq. (4.1) and decay branching 
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ratios in Fig. 3, are smaller than the measured branching ratios of corresponding non- 
resonant three-body decays of i?-mesons [28], specifically 

Br{B /C+7r°7r°) < 1.6 x 10~^ (4.2) 

Br(B K+TT^TT^ nonresonant) < 1.6 x 10"^, (4.3) 

BiiB ^ K+fi+fi-) kA.SxW'^, (4.4) 

Br (5 K+K^K+ nonresonant) < 2.8 x 10"^ (4.5) 

Resulting limits on model parameters are presented in Fig. 2. One can see, that a wide 
range of inflaton masses 300 MeV < < 5 GeV is allowed for direct investigation at B- 
factories. The branching ratio of the two-body i?-meson decay is at the level of 10~^ — 10"*^ 
for /3(^) in Fig. 1. 

Note in passing that our limits from i3-physics are only approximate. First, they 
do not apply to the inflaton masses close to the masses of known hadronic resonances, 
since limits (4.2)-(4.5) are obtained from the analysis of corresponding Dalitz plots with 
removed resonances. Second, the absence of new resonances on the Dalitz plots (would it 
be specially studied) put stronger constraints on the light inflaton model, than what we 
naively extracted from (4.2)-(4.5). Third, one should be careful about the bounds from 
B — )• Kfifi, because it is not valid for living inflatons > 10~^, which can escape the 
detector (the lower region near the CHARM bound). 



5. Discussion 



To summarize, following the new experimental results both in particle physics [8, 9] and 
cosmology [10, 12, 11], we have extended our previous analysis [1] of the large-field chaotic 
inflationary model with quartic inflaton potential to the case of inflaton non-minimally 
coupled to gravity. We concluded, that the range of cosmologically acceptable masses for 
the inflaton is wider, while the coupling between the inflaton and the SM sector can be 
stronger, than suggested [1] for the minimally coupled case. The model can be further tested 
with new cosmological data, especially with further results on the tensor perturbations of 
the CMB. 

The best place to probe the model remains the same: studies of two-body S-meson 
decays to kaons and inflatons, where the inflaton decays further into electron, muon, pion 
or kaon pair. However with the increase of the non- minimal coupling ^ the inflaton lifetime 
decreases. In a smaller part of parameter space the inflaton lifetime is long enough to cover 
a macroscopic distance from the i?-meson decay position, which leads to the following 
signatures of the inflaton event: either missing energy (if the inflaton escapes from the 
detector) or significantly displaced vertex of the inflaton decay. Therefore, searches for 
process B ^ K + nothing with two-body kinematics and B ^ K + x with subsequent 
decay of x to a pair of SM light particles at a macroscopic distance from the S-meson 
decay are necessary to probe the model. In a wider region of parameter space the inflaton 
decays close to the i?-meson decay point, thus effectively contributing to the three-body 
B-meson decays. This gives a possibility to probe the model by investigating neutral two- 
body resonance-like events {fi^fi~, vr+vr", etc.) among three-body final states of 5-meson 
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decays. The typical distance to the inflaton decay vertex, distinguishing these two options, 
can be deduced from the lifetime of the inflaton, see Figs. 2, 3. 

Independently (and also if S-meson decays are kinematically forbidden) the model 
can be probed to some extent in a fixed-target experiment. The inflaton can be produced 
by protons on target either in meson decays or through the gluon fusion similar to the 
SM Higgs boson, would it be of the same mass. This option has been studied in detail 
in Ref. [1] for the case of the inflaton minimally-coupled to gravity. The estimates of the 
inflaton production cross sections for a realistic set of beams are given there in Fig. 4. 
For the non-minimally coupled inflaton those estimates have to be multiplied by a factor 
/3(^)/(1.5 X 10^^^), which scales from 1 to 3 x 10^ when grows from to 1, see Fig. 1. 

Note also, that the model depends on the positivity of the SM Higgs boson self-coupling 
for the whole range of energy scales, relevant for inflation (or, to be more precise, for the 
scales of the order ~ ^ ajXMp). It implies a lower limit on the Higgs boson mass. This 
limit is slightly weaker than the bounds from the positivity of the Higgs self-coupling up to 
the Planck scale (see [29, 30, 31]), and is compatible with the observations at LHC within 
presently achieved precision for the Higgs and top-quark masses. 

Finally, the inflationary model under discussion can be further modified (completed) 
to solve other phenomenological problems of the SM (neutrino oscillations, dark matter 
phenomena, baryon asymmetry of the Universe) in such a way that the inflationary dy- 
namics and the low-energy phenomenology of the light inflaton remains unchanged. The 
working example is an extension of the SM with three Majorana sterile neutrinos, which 
can explain neutrino oscillations via type-I seesaw mechanism, the baryon asymmetry via 
resonant leptogenesis in the early Universe [32, 14] and the lightest of the sterile neutrinos 
can serve as dark matter. With sterile neutrinos in GeV mass range the model is known as 
i/MSM [13, 14] and thoroughly studied [33], though both seesaw mechanism and resonant 
leptogenesis allow for (much) heavier neutrinos as well. Note, that within our logic of only 
one energy scale in the whole model advertised in Introduction and Sec. 2, the natural 
mass scale for these neutrinos should not exceed the Electroweak scale. 

Moreover, the inflaton being singlet with respect to the SM gauge group may have 
Yukawa coupling to the sterile neutrino and hence contributes to their production in the 
early Universe and their masses after getting non-zero vacuum expectation value. In this 
way one can avoid the explicit mass term for the sterile neutrino and deal with the explicitly 
one-scale model of particle physics. This particular setup has been considered in Appendix 
A of Ref. [1] and all formulas presented there can be fully transferred to the case of the 
non-minimally coupled to gravity inflaton, since all of them explicitly contain the factor 
/3 = /3(^) and the inflaton mass. Most importantly, in the non-minimally coupled case 
as well, the inflaton may contribute to the dark matter production (the lightest sterile 
neutrino) thus eluding the strong flne-tuning [34, 35] required in the simplest version of 
i/MSM to generate the dark matter. Remarkably, the Yukawa interactions of inflaton with 
sterile neutrinos is bounded from the requirement of small radiative corrections to the 
inflationary potential, and can not give rise to very large masses [1, 7] consistently with 
the logic adopted. 

However, explicit mass parameters in the sterile neutrino part of the action are gener- 



-14- 



ally allowed. The only theoretical bound on them comes from the smallness of generated by 
the seesaw couplings one-loop finite corrections to the SM Higgs boson mass. The bound 
allows the sterile neutrino masses to be well above the EW scale. In this setup the mecha- 
nism of generation of the dark matter sterile neutrino in the inflaton decays still works, as 
well as the resonant leptogenesis [32]. We do not study this setup further as one beyond 
our logical framework. 
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